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Abstract 

We report current progress on the synthesis of methods to alleviate two ma- 
jor difficulties in implementing a Monte Carlo Renormalization Group (MCRG) 
for quantum systems. In particular, we have utilized the loop-algorithm to re- 
duce critical slowing down, and we have implemented an MCRG method in 
which the symmetries of the classical equivalent model need not be fully under- 
stood, since the Renormalization Group is given by the Monte Carlo simulation. 
We report preliminary results obtained when the resulting MCRG method is 
applied to the d=l XXZ model. Our results are encouraging. However, since 
this model has a Kosterlitz-Thouless transition, it does not yet provide a full 
test of our MCRG method. 



1 Introduction 

Although the Monte Carlo method introduced by Metropolis et al.^ has been a 
useful tool in the study of the critical behavior of classical systems, its application to 
the study of the critical behavior of quantum systems requires that many difficulties be 
circumvented. Some of these difficulties are already present in the study of classical 
systems, while others are present only in the study of quantum systems. In this 
paper we utilize present knowledge in overcoming the difficulties in such studies, and 
we apply the method to the study of the one- dimensional XXZ quantum spin model. 

* To appear in Quantum Monte Carlo Methods in Condensed Matter Physics, ed. M. Suzuki, 
World Scientific, 1993 
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The first difficulty to be surmounted in the study of quantum critical behavior 
was the problem of how to apply the Monte Carlo method to quantum systems. 
This difficulty was circumvented when Suzuki 0, |3|, |J] proposed to use the formula of 
Trotter || to map any rf-dimensional quantum system onto a <i+l-dimensional classical 
equivalent. Convergence properties of the Suzuki- Trotter transformation^) and of 
higher-order decompositions have been described Although this decomposition is 
applicable in general, it can lead to other difficulties not encountered in the study of 
most classical systems. 

One of the problems often introduced by the application of the Suzuki- Trotter 
transformation is the minus-sign problem. The minus-sign problem occurs when the 
'Boltzmann weights' of the classical equivalent of a quantum problem are not all pos- 
itive. In the Monte Carlo method Boltzmann weights are interpreted as probabilities, 
and the minus-sign problem to date is a serious one and can only be decreased. This 
is done by performing a simulation with any chosen probability distribution, and 
reweighting to the Boltzmann weights of the classical equivalent [§, ^ [H], [H], [T2|, [13 



Unfortunately, this moves the difficulty from one of principle to one of acquiring ad- 
equate statistics. In general the desired quantities are ratios of two numbers, each of 
which come from the difference between two large numbers with statistical errors. The 
minus-sign problem is therefore still largely unsolved. It is interesting to note that 
it is possible to make some classical statistical mechanical models have a minus-sign 
problem. For example, the simulation of the d—2 Ising ferromagnet has a minus-sign 
problem in a certain representation || |l4fl . 

A difficulty which occurs in the study of both classical and quantum critical be- 
havior is the critical slowing down that occurs in the standard Monte Carlo method 
near the critical point. The difficulty is that standard Monte Carlo algorithms em- 
ploy a local update procedure, and consequently 'information' diffuses through the 
lattice slowly in a random walk fashion. This difficulty was first mastered for the fer- 
romagnetic Potts model by Swendsen and Wang|l5| where the important realization 
was to perform non-local updates on clusters closely related to the critical clusters of 
the model. Similar procedures have been devised for a number of systems and have 
recently been reviewed fl6| , |T7f . The classical equivalent obtained from the Suzuki- 
Trotter transformation often maps onto a vertex model, and it is only recently that 
algorithms to alleviate critical slowing down for vertex models have been devised 

© HE H3 • 

In the study of phase transitions universal quantities, such as critical exponents, 
are desired. However, a phase transition occurs only in the thermodynamic limit - 
whereas Monte Carlo simulations by necessity can be done only on finite lattices. The 
best ways to obtain estimates of the critical exponents from Monte Carlo studies are 
to use finite-size scaling methods (for reviews see [23] , |22|[ ) or renormalization group 
(RG) methods|2^, p4| . One marriage of Monte Carlo and RG methods was given 
by Swendsen p5|, |26| , |2"T| . In Swendsen's implementation of the Monte Carlo Renor- 
malization Group (MCRG) one calculates only correlation functions (not coupling 
constants). To successfully apply this MCRG method, it is important that the RG 
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chosen preserves the symmetries of the Hamiltonian [28]. However, this presents a 
problem for quantum MCRG calculations since typically after the application of the 
Suzuki- Trotter formula one does not generally know the underlying 'Hamiltonian' of 
the d+l-dimensional classical equivalent. If you do not know the 'Hamiltonian', how 
can you know its symmetries? In particular, if the minus-sign problem is present, 
the concept of a 'Hamiltonian' on the classical equivalent is not well defined. If the 
classical equivalent is a vertex model, the constraint due to the continuity of vertex 
loops is difficult to preserve in an RG treatment. This is the difficulty that has de- 
layed progress on the application of MCRG procedures to quantum systems after the 



initial success of quantum MCRG methods to the transverse Ising model[29, 30, 31 



for which the underlying classical equivalent Hamiltonian is known and has no global 
constraints. For a review of quantum MCRG see Sec. 8 of Ref. [j3"2 ]. Similar diffi- 



culties have been encountered in applying Wilson's idea of real-space blocking p4| to 
exact-diagonalization studies, which has only recently been surmounted |3"B|, |3~3| , |55| . 

A method to overcome the difficulty of deciding what RG procedure to use in 
quantum studies was recently reported by Munger and Novotny|l4|. The idea is 



to modify a momentum-space MCRG developed by Swendsen[^7[], but to let the 
Monte Carlo simulation itself decide which RG to use. In momentum space this 
leads to the RG procedure of systematically discarding the least important degrees 
of freedom as one renormalizes the lattice. This momentum-space MCRG procedure 
was successfully applied to the square lattice g-state Potts model with q not restricted 



to integer values. Thus the MCRG in Ref. [|T4]] was performed on a staggered 6- vertex 
model which is obtained from a mapping of the g-state Potts model. 

In this paper we apply this MCRG for the first time to a quantum system — the 
6-vertex model, which is the classical equivalent of the d=l spin-| quantum XXZ 
model. The XXZ Hamiltonian is 



TV 



(i) 



with partition function Z = exp (—f37i). Periodic boundary conditions on the chain of 
N spins are used. Notice that 7i is invariant under (/3 — > —(3, A — > —A) by rotating 
every second spin through an angle of tt, and is thus equivalent to the Hamiltonian 



. crtcr. 



i+1 



+ a?a? +1 - Xa?af +1 



The case A=l corresponds to the isotropic Heisenberg antiferromagnet, while the 
case A=0 corresponds to the quantum XY chain. 

We have chosen this model since both the loop-algorithm of Ref. |T8| , [191 , pO|j , and 
the MCRG program of Ref. |14| were written to study the 6-vertex model, which is 



the classical equivalent of the Hamiltonian in Eq. ([!]) . 
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(See g%/, fig. 8.2) 



Figure 1: The 6 vertices which are allowed in the model are shown. For the XXZ model the weights 
are w$=ws=l and wi=W2=w^=w^. 



2 Monte Carlo Method 



It is well known ||, |37j, that the d—1 XXZ chain can be mapped to a 
(1+1) dimensional classical spin system by Suzuki's method. The main steps are a 
breakup of the Hamiltonian for an iV-spin chain into pieces living on even and odd 
bonds, an application of the Trotter-Suzuki formula, and an insertion of complete 
sets of intermediate states. The result is a system of classical spins Sij = ±1 on a 
periodic N x M chessboard lattice with the Euclidean time (Trotter) direction having 
extension M, and four-spin interactions on all "black" plaquettes of the chessboard 
[0, [], [SI], The mapping becomes exact as M — ► oo. 

The four spins (sjj; Sj+ij; Sij+i! Sj+ij+i) bordering an interacting plaquette obey 
the continuity constraint 



H,j+1 



+ S 



i+lj+l- 



(2) 



The world lines connecting sites with Sij = 1 are thus continuous. 

This classical spin system is a vertex model: vertices are located at the center of 
interacting plaquettes; they are connected by lines. Each such line touches a lattice 
site If Si j =+1 (—1), then we place an arrow that points upwards (downwards) 

in the Trotter direction on this line. Note that the resulting lattice of arrows and 
vertices has lines tilted 45 degrees w.r.t. the chessboard lattice. 

The constraint of Eq. (|2|) is that of the 6-vertex model: each vertex has two 
arrows pointing into the vertex and two pointing away from the vertex. The six 
possible vertex configurations are shown in Fig. 1. Their weights Wi in the partition 
function are the weights of the interacting plaquettes, namely |37[ 



a = W\ 
b = w 3 
c = w 5 



w 2 
w 6 



-A/3 



e A/9 sinh2/3, 
e A/9 cosh 2/3, 



(3) 



with (3 = [3/(2M). The 6-vertex model is exactly solved [42|. It is governed by the 
parameter 



A = 



a 2 + b 2 - c 2 



sinh(-2A/3) 3^0 



-A. 



(4) 



2ab sinh2/3 

Here we use periodic boundary conditions for the vertex lattice. The normal XXZ 
boundary conditions would be periodic in both the spin chain direction and the 



MCRG for d=l XXZ Model 



5 



Trotter direction, and correspond to periodic boundary conditions in the 45-degree 
directions for the vertex lattice. 

Monte Carlo simulations of constrained systems like the 6- vertex model have suf- 
fered from severe critical slowing down in the past. The situation has been drastically 
improved by the advent of the loop cluster algorithm [ |TB| , |T5| , |2"U| , which we now sketch 
very briefly. In terms of arrows, the constraint of Eq. (|2]) is a "divergence = zero" 
condition. All arrows therefore lie on directed closed loops (like magnetic field lines). 
The loop algorithm stochastically constructs such a loop and flips the direction of 
all arrows on the loop in a single Monte Carlo step while maintaining detailed bal- 
ance. It can therefore perform large steps in phase space, and produce statistically 
independent configurations within a few Monte Carlo steps even at infinite correla- 
tion length on large lattices [fT8| . Note that the loop algorithm also easily produces 
loops that wind around the lattice. At the KT transition (- = ~ = |), for example, 
autocorrelation times are reduced by about two orders of magnitude on a 64 2 lattice 
with respect to a Metropolis simulation. 



3 MCRG Method 

Here we use the MCRG method described in detail in Ref. |14|]- The idea of this 
MCRG method is to let the Monte Carlo simulation itself provide the renormaliza- 
tion group to be used. This is done by building on the MCRG method developed by 



Swendsen [p5| , p26fl to calculate critical exponents. Rather than use a block spin pro- 
cedure in real space, the RG is defined in momentum space in a similar fashion as the 
procedure introduced by Swendsen in Ref. pTI] . Rather than using the Fourier trans- 



form of the configurations [27], we use the Fourier transform of each state individually 

#k = H <Wr exp (ik • r) , (5) 

r 

where a labels the vertex state v, and the summation is over the simulated lattice. In 
our case, a takes the values of one through six. An inverse Fourier transform over a 
restricted part of the momentum-space is then performed. The same restricted part 
is used for all generated configurations, and the inverse transform is given by 

k 



C = E'^exp -i(k-m Q ). S , (6) 



where is a shift in momentum space, and the prime on the summation indicates 
that the sum is over a restricted part of momentum space. Finally, on each site of the 
reduced lattice for each configuration the state a is chosen for which the real part, 
9ft(t>g a ) is largest. Of course this choice is arbitrary, since for example the state with 
the largest modulus could have been chosen. 

Both the restricted part of momentum space and the shift m a are the same for 
all generated configurations and are given by the Monte Carlo simulation itself. The 
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shift rrio, is chosen such that the important fluctuations, as determined by the peak 
in the momentum-space plot from the entire Monte Carlo simulation, are shifted to 
the center of the Brillouin zone. This transforms each state into a mainly 'ferromag- 
netic' state, and allows the use of normal ferromagnetic operators to obtain critical 
exponents from the linearized transformation matrix T*. The restricted region of the 
shifted momentum space is chosen to be the region in the center of the Brillouin zone 
that gives N'—L'xL' spins on transforming back to real space. Note that this MCRG 
procedure does not preserve the constraints of the vertex model (the continuity of 
loops). However, the renormalized variables should be considered to be composite 
variables, and the vertex constraints on the original lattice are taken into account 
during the RG procedure in a way dictated by the Monte Carlo simulation of the 
model. 

In principle, one could avoid transforming back to real space, but at the expense 
of defining the 'majority rule' for the states in momentum space, and at the expense 
of defining operators for the linearized transformation matrix T* in momentum space; 
neither of which has been attempted to the authors' knowledge. The transformation 
back to real space reduces the lattice size and gives the RG. However, a normalization 
is also required. This is done by multiplying the new real-space spins with a constant 
such that J2\ v r\ 2 /N = Yl'\ v ' s \ 2 /N', where N' is the number of lattice sites on the 
reduced lattice, and the primed summation is over renormalized spins v' s on the 
reduced lattice. Other normalization methods can be devised |TJ]], but they should 
all lead to the same critical exponents from the MCRG. 

The critical properties are obtained in the usual fashion from the linearized RG 
transformation matrix T* from the elements 

dK (n) 

T* — — — 3 (7) 

a/3 a T A n + m ) ' 1 ' 

where is the coupling constant corresponding to the operator S a following the 
n th RG transformation. The K a include factors of 1/ksT. From the chain rule, 

d(s^) d(sy>) 

where the angular brackets indicate the thermal average of the operator. The partial 
derivatives can be easily calculated from the canonical ensemble: 

= {S ^st +m) ) - (S^)(S£ +m ~>). (9) 

OKa 

The thermal averages in Eq. (Q) are calculated between the renormalized lattices after 
n and n+m RG transformations. These averages can be calculated at the simulated 
set of parameters, but they can also be calculated with parameters other than the 
simulated ones using a single reweighting method f|, |K], |TT| |12| , or using data from 



multiple simulations for different parameters [12, 13, 14 . 
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(See gjy, fig. 8.5) 



Figure 2: The phase diagram for the 6-vertex model is shown. The vertex weights satisfy W1—W2, 
and u>5=u>6- The dashed line corresponds to the F-model where W1—W3 as well, and this 
line terminates at the 'infinite temperature' point where W1—W5 (filled circle). The open circles 
represent the points where we have made MCRG calculations. 



The number of operators in Eq. (R) is infinite. However, the MCRG estimate for 
the matrix T* is found from Eq. (§) by first truncating all the matrices in Eq. (§) and 
then multiplying Eq. (||) by the inverse of the truncated matrix that multiplies T* in 
Eq. (|8]). To first order in perturbation theory, it has been shown that this procedure 
takes into account the operators which have not been included in the truncation 
procedure |39[]. Since we have performed the shift m a for each state, we need to 



include only 'ferromagnetic' operators. 

Once the largest eigenvalues, Ai>A2>- • •, of the linearized transformation matrix 
T*^ are determined, the eigenvalues of the RG are given by 

y 4 = ln(A,)/ln(6), (10) 

where the scale factor b is given by the ratio of the lengths on the two lattice sizes used, 
b=L f iargJ ' L 'smaU={ N L g J Kmaii) 1/d > where d is the dimension of the model studied. 
The exponents yi are irrelevant if they are less than zero, are marginal if Vi=0, and 
are relevant and give universal critical exponents for models associated with the fixed 
point if yi>0. 

For a second-order transition the two largest critical exponents of the RG (yr from 
the operators in the thermal section and yn from operators in the magnetic section) 
allow the critical exponents to be calculated in the normal fashion from v=l/yT and 
r]=2+d-2y H . 

Special care must be taken in the case of a model with a Kosterlitz-Thouless (KT) 



transitional . The authors in Ref. [[5J show that although the magnetic exponents 
are given by r)=2+d—2yH in the normal fashion, the values of v for a KT transition 
must be obtained by fitting to the KT RG equations. We have not done such a 
fitting procedure, but have rather concentrated only on the magnetic exponents. 
In general one can also use the MCRG method to obtain the values of the critical 
couplings |40"|, f41~l . However since the critical couplings for the d=l XXZ model are 
known exactly, we have not tried to find them using MCRG. 



4 Theoretical results for the d=l XXZ Model 



The d=l spin-i XYZ model has been studied using a wide variety of meth- 
ods. These methods have ranged from exact methods[52, 44], to universality 
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and mapping methods [45, 46], to conformal invariance methods [47, 48 



to exact diagonalization studies [|50 



to series ex- 



52 , to Monte Carlo 



pansion methods f49 

methods 0, |5^, [54]. The phase diagram which has emerged from these studies is il- 
lustrated in the 6- vertex representation in Fig. 2 [ 42 . We have performed our MCRG 
calculations only on the F-model (where Wi=W2=w 3 =w 4: and w 5 =w 6 ) at points illus- 
trated by the open circles in Fig. 2. In Fig. 2 region I corresponds to a ferroelectric 
region where the lowest energy state is one with all vertices equal to 1 or 2. In this 
region the excited states give a negligible contribution to the partition function and 
the system is frozen in one of the two ground states. The situation is similar for 
region II, where all vertices equal 3 or 4, and region IV where the two ordered states 
have vertices of 5 on one sublattice and 6 on the other sublattice. Region III is the 
disordered phase where there is no spontaneous order, nor interfacial tension, but the 
correlation length is infinite and correlations decay as an inverse power of distance 
rather than exponentially [f4"2"|] . 

The F-model with wiKw^/2 (which corresponds to the asymmetry parameter of 
Eq. ([!]), A>1) is in region IV. Since in this region the model is completely ordered, 
one expects to obtain a low-temperature discontinuity fixed point (the Nienhuis- 
Nauenberg criterion |55|) with the largest RG eigenvalue y=d. 

There is a Kosterlitz-Thouless (KT) transition in the F-model when W\=w^j2 
(where A=l, and the model is the Heisenberg model). In the XXZ model there is 
a KT transition from the massless into the antiferroelectric ordered domain at A=l. 
This occurs through a Kosterlitz-Thouless mechanism caused by the excitations of an 
umklapp process[4~5", 16 . At A=l the vortex operator (which comes from the critical 
exponent of the 8-vertex model) is marginal. 

For the F-model in region III the critical exponents are those of the critical 8-vertex 
model|42|] and are the same as those of the Luttinger model[[J5[]. From conformal 
invariance ||48|| the entire operator content is known. Stringlike solutions of the Bethe- 
Ansatz equations yield excited states corresponding to operators Oj ,• with dimensions 



x 



i 2 x p + f/Axp i,j 



0,1,2, 



mi 



These operators are the analogs of the Gaussian-model operators composed of spin- 
wave excitations of index i and a "vortex" of vorticity j. In Eq. ([11]) x p =(tc— r y)/27T, 
where 7=arccos(A). The operators Oi,o and O2,o correspond respectively to the polar- 
ization and energy operators of the 8-vertex model. The operator Oo,i is the crossover 
operator of the 8-vertex model, or equivalently the energy operator of the Ashkin- 
Teller model. From the dimensions Xij the RG exponents are found through 

d — Xij. (12) 



The exponents i] i j=2+d—2y i j=2xij govern the behavior of correlations. For exam- 
ple, for o ,i one has 77o,i=2xo,i=2+<i— 2t/ ,i and this critical exponent governs the 
correlation 77,2=770,1 with 



{S z (0)S z (r)} ~ r 



(13) 



where S z =^a z 
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(See m], fig. 4) 



Figure 3: The 11 operators used for the truncated linearized RG transformation matrix T* are 
shown. Filled circles represent a vertex in any state, circles with numbers inside represent a S- 
function for that vertex state, and the 'bull's eye' in the last operator is any vertex which is not a 
5 or 6. The square lattice is drawn with light lines, while heavy lines join the interacting sites. For 
the operator to be nonzero all filled-circle interacting vertices must have the same vertex state. All 
combinations which are related by translation, rotation, and reflection to those shown were included 
in the calculation. 



5 MCRG results 

We have applied the MCRG method described in Sec. 2 to the 6- vertex model on 
LxL lattices with periodic boundary conditions. The 11 operators which we included 
in the truncated linearized RG transformation matrix T* are shown in Fig. 3. Note 
that because the momentum-space RG does not conserve the 6-vertex constraint on 
the renormalized lattice, we have not included polarization-type operators in the 11 
measured operators. Our calculations were done on a CRAY Y-MP432 supercom- 
puter, and the Fourier transforms used were FFT's in the NAG library (the explicit 
routines used depended on the lattice sizes the FFT was applied to). 

We first studied the F-model with in the ordered region IV of Fig. 2, 

and the results are shown in Table 1. The lattice simulated for Table 1 was 32x32 
and the renormalized lattices were L'xL' with L'=27, 25, 21, 17, and 15. The RG 
analysis was performed using 15000 configurations which were generated skipping 250 
Monte Carlo cluster updates between configurations and with 250000 cluster updates 
for thermalization. The exponential autocorrelation time, r exp in units of 'sweeps' |[20|| 
is given by the autocorrelation time in units of cluster updates times the average 
cluster size divided by 2L 2 . Here r exp =43(5), and it is most clearly visible in the 



sublattice energy p0|. The average cluster size is 0.96 times 2L 2 . In principle one can 
obtain estimates for the statistical errors of the eigenvalues in Table 1, for example 
by analyzing J bins of the generated configurations. We have not yet performed such 
an analysis. 

The values of Tables 1 and 2 should be read in a specific fashion to see the 
convergence of the exponents. In order to obtain good critical exponents from the 
MCRG one needs to be able to penetrate the linear region about the fixed point, since 
the exponents are calculated from the linearized transformation matrix T*^. It may 
take several iterations to be able to get close enough to the fixed point that a linear 
approximation is reasonable. Consequently, the first iterations (the ones starting from 
the largest lattices, and hence on the first lines in the tables) may need to be ignored, 
or at least to be used only to see how the convergence toward the linear region is 
proceeding. For this reason we do not show the iteration between the original lattice 
and the next largest lattice in the tables. The exponents should then converge for 
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1 L 


777 1 2 3 4 




J-Jdl^CoL 1\ VI Il/AUUIICIIL 


1 
1 




2 


2.0003 1.9980 1.9868 


3 


1.9963 1.9840 


4 


1.9811 




Second RG exponent 


1 


0.6449 1.0041 1.2214 1.1063 


2 


0.7962 1.1325 1.0128 


3 


0.9313 0.7794 


4 


0.9683 




RG scale factor b 


1 


1.080 1.286 1.588 1.800 


2 


1.191 1.471 1.667 


3 


1.235 1.400 


4 


1.133 



Table 1: The two largest exponents of the linearized transformation matrix, T*, are listed for the 
simulation in the ordered phase at wi/w§=\. The lattice simulated was a 32x32 lattice, and the 
sizes of the renormalized L'xL' lattices were L'—27, 25, 21, 17, and 15. Listed are the exponents 
between renormalized lattices n and n+m as well as the corresponding RG scale factors b. For 
example, the entry with n=2 and m=3 corresponds to an RG transformation between lattices with 
L'=25 and L'=15. The data should converge for large n,m, but also become susceptible to finite 
size effects. The expected largest eigenvalue is that of the zero-temperature fixed point since the 
system is ordered in Region IV of Figure 2. Consequently the largest exponent should be that of 
the discontinuity fixed point, y=d=2. Our data converge to this value extremely rapidly. The next 
largest eigenvalue would be the first analytic correction, which approaches unity. 



large values of n and m. 

If n, m become too large, however, there are additional difficulties. Smaller scale 



factors b will have the larger statistical errors[|T4|. Also, in MCRG studies the larger 



statistical errors typically occur on the smallest lattices, which correspond to the last 
entry in each line of the tables. The smallest lattices will also have larger finite-size 
effects, which are particularly disruptive for the longer-range operators included in 
the calculation of T*g. 

In Table 1 one sees extremely rapid convergence to the eigenvalue y=d=2 associ- 
ated with the Nienhuis-Nauenberg criterion for a discontinuity fixed point ||55|| . The 
next-largest eigenvalue, while not having converged too well, seems to be approaching 
the value of unity, which would give the first analytical correction to the discontinuity 
fixed-point behavior. 

Our other MCRG results are shown in Table 2. These are at the KT point 
(A=l), the XY point (free-fermion case with A=l/ v^2), and an intermediate point with 
A=l/2. Table 2 lists the two largest eigenvalues obtained from our MCRG procedure. 
Each value of A studied took approximately 23 hours of Y-MP time for the RG analysis 
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Wl /w 5 = l/2 A=l 




Wl /w 5 = l/y/2 A = 


n 


m 1 2 3 4 


n 


m 1 2 3 4 




Largest RG Exponent 




Largest RG Exponent 


1 


2.022 1.809 1.707 1.667 


1 


1.902 1.477 1.342 1.299 


2 


1.676 1.624 1.605 


2 


1.098 1.145 1.165 


3 


1.591 1.581 


3 


0.975 1.070 


4 


1.563 


4 


0.844 




Second RG Exponent 




Second RG Exponent 


1 


1.466 1.068 0.835 0.763 


1 


1.185 0.493 0.175 0.198 


2 


0.636 0.585 0.576 


2 


-.520 -.234 -.005 


3 


0.436 0.500 


3 


-.610 -.166 


4 


0.358 


4 


-.931 





Wi/iyg = 0.5825 •• • A = 1/2 


n 


m 1 2 3 4 




Largest RG Exponent 


1 


2.050 1.715 1.581 1.524 


2 


1.501 1.463 1.437 


3 


1.397 1.391 


4 


1.305 




Second RG Exponent 


1 


1.429 0.908 0.616 0.537 


2 


0.292 0.307 0.327 


3 


0.029 0.216 


4 


-.137 


n 


m 1 2 3 4 




RG scale factor b 


1 


1.244 1.645 2.429 3.400 


2 


1.323 1.952 2.733 


3 


1.476 2.067 


4 


1.400 



Table 2: The two largest exponents of the linearized transformation matrix T* are listed for the 
simulations at A=0, \ 1 and 1. The lattice simulated was 64x64, and the sizes of the renormalized 
L'xL' lattices were Z/=51, 41, 31, 21, and 15. Listed are the exponents between renormalized 
lattices, n and n+m, as well as the corresponding RG scale factors b. For example, the entry with 
n=3 and m=2 corresponds to an RG transformation between lattices with L'=31 and L'=15. The 
data should converge for large n, m; but should also become susceptible to finite size effects. For 
A=0, \, and 1, the largest eigenvalues y z =2—r] z , corresponding to the operator Oo,i, are 1, | and 
|, respectively. 
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of the generated configurations. The time required for the Monte Carlo portion of 
the runs depended on the value of A since the cluster-size per volume depended on 
A. For A=0, |, and 1 the cluster-size was 0.15, 0.09, 0.038 times 2L 2 and the Y-MP 
CPU time used corresponded to 3.9, 2.5, and 1.4 hours. The lattice simulated for 
Table 2 was 64x64, and the renormalized lattices were L'xL' with L'=51, 41, 31, 
21, and 15. The RG analysis was again performed using 15000 configurations which 
were generated skipping 250 cluster updates between configurations and using 250000 
cluster updates for thermalization. Consequently, a total of 4xl0 6 cluster updates 
were generated. The exponential autocorrelation time, which is most clearly visible 
in the sublattice energy |2"0||, for A=0, |, and 1, is r exp =4.0(2), 8.1(9), and 5.7(6), 
respectively. 

The convergence toward the linear region about the fixed point is expected to be 
particularly slow when a marginal eigenvalue is present, as is the case in Table 2. 
Consequently, in Table 2 convergence has not yet occured, but the trend toward 
convergence of the largest eigenvalue is evident and in agreement with the expected 
result. 

In Table 2 we see only one nearly converged largest exponent y for each value of A. 
This exponent can be associated with the operator O ,i corresponding to the cross- 
over exponent of the 8- vertex model, and the exponent r] z can be associated with the 
S z correlations. The reason operators such as Oi o may not be found in our MCRG 
calculation is because such operators correspond to the polarization operator of the 
8-vertex model, and such asymmetric operators are not included in the operators 
shown in Fig. 3. Such operators are not easily included since the momentum-space 
RG procedure does not preserve the local vertex constraint. 



6 Discussion and Conclusions 

We have implemented for the first time the marriage of a loop algorithm to alle- 
viate critical slowing down and an MCRG procedure that does not necessitate that 
one know the symmetries of the underlying classical equivalent beforehand since they 
are given by the Monte Carlo. We have applied the resulting MCRG method to the 
study of the d—1 spin-| quantum XXZ chain. Although this model has a Kosterlitz- 
Thouless transition, and consequently the convergence toward the fixed point is fairly 
slow, we are able to obtain the critical exponent r] z associated with the S z correlation 
function reasonably well. The exponent r] z comes from the operator Oo,i, which is 
associated with the crossover exponent of the 8-vertex model. We have pointed out 
the difficulties in our study that have prevented us from obtaining the other expo- 
nents of the XXZ model. Although this is a preliminary study, our MCRG method 
shows a great deal of promise, since it is easily generalized to the study of many other 
quantum models in one and higher dimensions. 
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